Strong violation of critical phenomena universality: 
Wang-Landau study of the 2d Blume-Capel model under bond randomness 



A. Malakis^ A. Nihat Berker'^-^'^, I. A. Hadjiagapiou\ and N. G. Fytas^ 
^Department of Physics, Seetion of Solid State Physics, 
University of Athens, Panepistimiopolis, GR 15784 Zografos, Athens, Greece 
^ College of Sciences and Arts, Kog University, Sariyer 34450, Istanbul, Turkey 
^Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, U.S.A. and 
"^Feza Giirsey Research Institute, TUBITAK - Bosphorus University, Qengelkdy 34684, Istanbul, Turkey 

(Dated: March 1, 2009) 

We study the pure and random-bond versions of the square lattice ferromagnetic Blume-Capel 
model, in both the first-order and second-order phase transition regimes of the pure model. Phase 
transition temperatures, thermal and magnetic critical exponents are determined for lattice sizes 
in the range L = 20 — 100 via a sophisticated two-stage numerical strategy of entropic sampling 
in dominant energy subspaces, using mainly the Wang-Landau algorithm. The second-order phase 
transition, emerging under random bonds from the second-order regime of the pure model, has 
the same values of critical exponents as the 2d Ising universality class, with the effect of the bond 
disorder on the specific heat being well described by double-logarithmic corrections, our findings 
thus supporting the marginal irrelevance of quenched bond randomness. On the other hand, the 
second-order transition, emerging under bond randomness from the first-order regime of the pure 
model, has a distinctive universality class with v = 1.30(6) and P/v — 0.128(5). These results 
amount to a strong violation of universality principle of critical phenomena, since these two second- 
order transitions, with difi'erent sets of critical exponents, are between the same ferromagnetic and 
paramagnetic phases. Furthermore, the latter of these two sets of results supports an extensive but 
weak universality, since it has the same magnetic critical exponent (but a different thermal critical 
exponent) as a wide variety of two-dimensional systems with and without quenched disorder. In 
the conversion by bond randomness of the first-order transition of the pure system to second order, 
we detect, by introducing and evaluating connectivity spin densities, a microsegregation that also 
explains the increase we find in the phase transition temperature under bond randomness. 

PACS numbers: 75.10.Nr, 05.50.-|-q, 64.60.Cn, 75.10.Hk 



INTRODUCTION: STRONG VIOLATION OF 
UNIVERSALITY 



Universality, according to which the same critical ex- 
ponents occur in all second-order phase transitions be- 
tween the same two phases, erstwhile phenomenologically 
established, has been a leading principle of critical phe- 
nomena [l[ . The explanation of universality, in terms of 
diverse Hamiltonian flows to a single fixed point, has been 
one of the crowning achievements of renormalization- 
group theory [2|. In rather specialized models in spatial 
dimension d = 2, such as the eight- vertex Q and Ashkin- 
Teller 14] models, the critical exponents nevertheless vary 
continuously along a line of second-order transitions. We 
shall refer to these cases as the weak violation of univer- 
sality. We have established in the current study a much 
stronger and more general instance of universality viola- 
tion, under the effect of quenched bond randomness. It 
has been known that quenched bond randomness may or 
may not modify the critical exponents of second-order 
phase transitions, based on the Harris criterion Q. 
It was more recently established that quenched bond 
randomness always affects first-order phase transitions 
by conversion to second-order phase transitions, for in- 
finitesimal randomness in d — 2 [7, 8 and after a thresh- 
old amount of randomness in d > 2 [s | , as also inferred by 
general arguments [9| . These predictions 0, 3 have been 



confirmed by Monte Carlo simulations lid. More- 
over, renormalization-group calculations [12l | on tricriti- 
cal systems have revealed that not only first-order tran- 
sitions are converted to second-order transitions, but the 
latter are controlled by a distinctive strong-coupling fixed 
point. 

In the current Wang-Landau (WL) study yielding es- 
sentially exact information on the two-dimensional (2d) 
Blume-Capel model under quenched bond randomness, 
we find dramatically different critical behaviors of the 
second-order phase transitions emerging from the first- 
and second-order regimes of the pure model. These 
second-order transitions with the different critical ex- 
ponents are between the same two phases indicating a 
strong violation of universality, namely different sets of 
critical exponents on two segments of the same critical 
line. Moreover, the effect of quenched bond random- 
ness on the critical temperature is opposite in these two 
regimes, which we are able to explain in terms of a mi- 
crosegregation mechanism that we observe. Finally, in 
proving a general strong violation of universality under 
quenched bond randomness, our study supports a more 
delicate and extensive weak universality [ij, : In the 
random-bond second-order transition emerging from the 
pure-system first-order transition, the magnetic (but not 
thermal) critical exponent appears to be the same as that 
of the pure 2d Ising model, as has also been seen in other 
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random and non-random systems. 

The Blume-Capel (BC) model is defined by tlie 

Hamiltonian 

Hp^-JJ2^^SJ+^J2^^' (1) 
<ij> i 

where the spin variables St take on the values —1,0, 
or +1, < ij > indicates summation over all nearest- 
neighbor pairs of sites, and the ferromagnetic exchange 
interaction is taken as J = 1. The model given by Eq. ([T]), 
studied here in 2d on a square lattice, will be referred to 
as the pure model. Our main focus, on the other hand, 
is the case with bond disorder given by the bimodal dis- 
tribution 

P{J^l) = \ W^, - Ji) + - J2)] ; (2) 

J1+J2 . J 2 

= 1 ; r = — , 

2 ' Ji ' 

so that r reflects the strength of the bond randomness. 
The resulting quenched disordered (random-bond) ver- 
sion of the Hamiltonian defined in Eq. ^ reads now as 

<ij> i 

II. TWO-STAGE ENTROPIC SAMPLING 

We briefly describe our numerical approach used to es- 
timate the properties of a large number, 100, of bond 
disorder realizations, for lattice sizes L = 20 — 100. The 
pure-system properties are also obtained, for reference 
and contrast. We have used a two-stage strategy of a re- 
stricted entropic sampling, which is described in our re- 
cent study of random-bond Ising spin models in 2d p^j , 
very similar to the one applied also in our numerical 
approach to the 3d random-field Ising model [l^. In 
these papers, we have presented in detail the various 
sophisticated routes used for the restriction of the en- 
ergy subspace and the implementation of the WL al- 
gorithm 19, 20]. The identification of the appropriate 
energy subspace {Ei,E2) for the entropic samphng of 
each random-bond realization is carried out by applying 
our critical minimum energy subspace (CrMES) restric- 
tion [U [l^] and taking the union subspace at both pseu- 
docritical temperatures of the specific heat and suscepti- 
bility. This subspace, extended by 10% from both low- 
and high-energy sides, is sufhcient for an accurate esti- 
mation of all finite-size anomalies. Following Ref. p7j . 
the identification of the appropriate energy subspace is 
carried out in the first multi-range (multi-R) WL stage 
in a wide energy subspace. The WL refinement levels 
iG{E) f * G{E), where G{E) is the density of states 
(DOS); for more details see Ref. [13]) used in this stage 
U = l,---,ii;/i+i = i/ii) were ji = 18 for L < 80 
and ji = 19 for L > 80. The same process was re- 
peated several times, typically ~ 5 times, in the newly 



identified restricted energy subspace. From our experi- 
ence, this repeated application of the first multi-R WL 
approach greatly improves accuracy and then the result- 
ing accurate DOS is used for a final redefinition of the 
restricted subspace, in which the final entropic scheme 
(second stage) is applied. In this stage, the refinement 
WL levels j = ji, . . . , ji + A are used in a one-range (one- 
R) or in a multi-R fashion. For the present model, both 
approaches were tested and found to be sufficiently accu- 
rate, provided that the multi-R uses adequately large en- 
ergy subintervals. This fact will be illustrated in the fol- 
lowing section, by presenting the rather sensitive double- 
peak (dp) structure of the energy probability density 
function (PDF) in the first-order regime of the model. 
Noteworthy, that most of our simulations of the 2d BC 
model at the second-order regime (A = 1) were carried 
out by using in the final stage a one-R approach, in which 
the WL modification factor was adjusted according to 
the rule In / ~ proposed recently by Belardinelli and 
Pereyra [23] . Our comparative tests showed that this al- 
ternative approach yields results in agreement with the 
one-R WL approach. 

Let us close this brief outline of our numerical scheme 
with some appropriate comments concerning statistical 
errors and disorder averaging. Even for the larger lat- 
tice size studied here (L — 100), and depending on the 
thermodynamic parameter, the statistical errors of the 
WL method were found to be of reasonable magnitude 
and in some cases to be of the order of the symbol sizes, 
or even smaller. This is true for both the pure version 
and the individual random-bond realizations. These WL 
errors have been used for the pure system in our finite- 
size scaling (ESS) illustrations and fitting attempts. For 
the disordered version only the averages over the disor- 
der realizations, denoted as [. . .]av, have been used in the 
text and their finite-size anomahes, denoted as [. . .]*^,, 
have been used in our ESS attempts. Due to very large 
sample-to-sample fluctuations, mean values of individ- 
ual maxima ([. . .*]av) have not been used in this study. 
However, for the finite-size anomalies of the disordered 
cases, the relevant statistical errors are due to the finite 
number of simulated realizations. These errors were es- 
timated empirically, from runs of 20 realizations via a 
jackknife method, and used in the corresponding ESS fit- 
ting attempts. These disorder-sampling errors may vary, 
depending again on the thermodynamic parameter, but 
nevertheless were also found to be of the order of the 
symbol sizes. For the case A = 1, these are hardly ob- 
servable, as illustrated in the corresponding graphs. 



III. PHASE TRANSITIONS OF THE PURE 2d 
BC MODEL 

A. First-order transition of the pure model 

The value of the crystal field at the tricritical point of 
the pure 2d BC model has been accurately estimated to 
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FIG. 1: (color online) Behavior of the pure 2d BC model at A = 1.975: (a) FSS behavior of the specific heat peaks giving 
a clear divergence, (b) The same for the susceptibility maxima, (c) Simultaneous fitting of the maxima of the averaged 
logarithmic derivatives of the order parameter defined in Eq. (j4j. (d) Power- law behavior of the averaged absolute order- 
parameter derivative in Eq. ([5]). (e) The dp structure of the energy PDF at T = Th via the two difi'erent implementations 
of the WL scheme. The multi-R implementation is displaced very slightly to the right, (f) Limiting behavior for the surface 
tension E(L) defined in the text and the latent heat Ae(L) shown in panel (e). The linear fits shown include only the five 
points corresponding to the larger lattice sizes. 



be At = 1.965(5) [H, [H, [H, [H . Therefore, we now 
consider the value A = 1.975, for which the pure model 
undergoes a first-order transition between the ferromag- 
netic and paramagnetic phases, and carry out a detailed 
FSS analysis of the pure model. Our first attempt to elu- 
cidate the first-order features of the present model will 
closely follow previous analogous studies carried out on 
the g = 5, 8, 10 Potts model [29|, M, M, US, HI and also 
our study of the triangular Ising model with competing 
nearest- and next-nearest-neighbor antiferromagnetic in- 
teractions [S^]. As it is well known from the existing 
theories of first-order transitions, all finite-size contribu- 
tions enter in the scaling equations in powers of the sys- 
tem size L"^ [IHl . This holds for the general shift behavior 
(for various pseudotransition temperatures) and also for 
the FSS behavior of the peaks of various energy cumu- 



lants and of the magnetic susceptibility. It is also well 
known that the dp structure of the energy PDF, P{e), 
where e = H/L^, is signaling the emergence of the ex- 
pected two delta-peak behavior in the thermodynamic 
limit, for a genuine first-order phase transition [sgI. \3T\. 
and with increasing lattice size the barrier between the 
two peaks should steadily increase. According to the 
arguments of Lee and Kosterlitz [s^, HH the quantity 
AF{L)/L = [kBTln{P^ax/Prmn)]/L, where P„,ax and 
Pmin arc the maximum and minimum energy PDF val- 
ues at the temperature Th where the two peaks are of 
equal height, should tend to a non-zero value. Similarly 
to the above, the logarithmic derivatives of the powers of 
the order parameter with respect to the inverse temper- 
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FIG. 2: Behavior of the pure 2d BC model at A = 1: (a) Simultaneous fitting of six pseudocritical temperatures defined in 
the text for L > 50. (b) FSS of the specific heat peaks. The inset shows a linear fit of the specific heat data on a log scale 
for L > 50. (c) FSS of the susceptibility peaks on a log-log scale, (d) FSS of the order parameter at the estimated critical 
temperature also on a log-log scale. Linear fits are applied in panels (c) and (d). 
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and the average absolute order-parameter derivative, 
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have maxima that scale as L"^ with the system size in the 
FSS analysis of a first-order transition. In the case of 
a second-order transition, the quantities in Eqs . (ffl and 
© respectively scale as L^/"" and L'^^-f^)/'' [13, O HI] , 
to be used further below. 

Figures [IJa) and (b) illustrate that the traditionally 
used divergences in FSS of the specific heat C and suscep- 
tibility X follow very well a power law of the form ^ L"^, as 
expected for first-order transitions [s^, ■ Furthermore, 
Figs. [T{c) and (d) demonstrate that the divergences cor- 
responding to the first-, second-, and fourth-order log- 
arithmic derivatives of the order parameter defined in 
Eq. dH) and the absolute order-parameter derivative de- 
fined in Eq. (O follow also very well the same i'^ behav- 
ior, as expected. Figure [T]^e) shows the pronounced dp 
structure of the energy PDF of the model at T = Th for 
L = 60, obtained by the two different implementations of 
the WL scheme. This graph illustrates that the dp struc- 
ture is not very sensitive to the multi-R WL process, in 
contrast to our recent findings for some first-order-like 
behavior of the 3d random-field Ising model. It also il- 
lustrates the accuracy of the implementation schemes. 
As mentioned above, from these dp energy PDF's one 
can estimate the surface tension = AF{L)/L and 

the latent heat Ae(L), whose values remain finite for a 
genuine first-order transition. Figure [Iff) shows the lim- 
iting behavior of these two quantities and verifies the 
persistence of the first-order character of the transition 
at A = 1.975. The limiting values of S(L) and Ae(L) 
are given in the graph by extrapolating at the larger lat- 
tice sizes studied. We close this section by noting that 



the transition temperature T*(A = 1.975) is estimated 
to be, in the limit L ^ oo, T* = 0.574(2). This value 
interpolates and agrees with the general phase diagram 
points summarized in Ref. [28j . 



B. Second-order transition of the pure model 

The 2d BC model with no quenched randomness, 
Eq. at the crystal field value A = 1, undergoes a 
second-order transition between the ferromagnetic and 
paramagnetic phases, expected to be in the universality 
class of the simple 2d Ising model [1^. In the follow- 
ing, we present the FSS analysis of our numerical data 
for this case, to verify this expectation and to set any 
contrast with the behavior under quenched randomness, 
presented further below. Figure [2lja) gives the shift be- 
havior of the pseudocritical temperatures corresponding 
to the peaks of the following six quantities: specific heat, 
magnetic susceptibility, derivative of the absolute order 
parameter, and logarithmic derivatives of the first, sec- 
ond, and fourth powers of the order parameter. Fitting 
our data for the larger lattice sizes (L = 50 — 100) to 
the expected power-law behavior T — Tc + bL~^/'^ , we 
find that the critical temperature is Tc — 1.3983(5) and 
the shift exponent is — 1.013(14). Almost the same 
estimate for the critical temperature is obtained when 
we fix the shift exponent to the value 1/v = 1. Thus, 
the shift behavior of the pseudocritical temperatures in- 
dicates that the pure 2d BC model with A = 1 shares 
the same correlation length exponent v with the 2d Ising 
model. Figure [U^b) gives the FSS of the specific heat 
peaks. Here, the expected logarithmic divergence of the 
specific heat is very well obtained even from the smaller 
lattice sizes as shown in the main frame. The inset is 
a linear fit of the specific heat data on a log scale for 
L > 50. Finally, Figs, [^c) and (d) present our estima- 
tions for the magnetic exponent ratios "//v and f3/v. In 
panel (c) we show the FSS behavior of the susceptibility 
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FIG. 3: (color online) Behavior of the random-bond 2d BC model at A = 1.975: (a)-(d): Softening effects on the first-order 
transition features of the 2d BC model, induced by bond randomness of various strengths. In (a) and (b), the very rough P(e) 
is for the pure model, whereas the smoothed curves are for r — 0.90 (deeper) and r = 0.82 (shallower). The curves in (c) and 
(d) are for different seeds. Panels (e) - (h) illustrate the (e) - T behavior and (M) - T behavior, for the size L = 30 and various 
disorder strengths r for different disorder realizations. The curves in (e) and (f) are for r = 0.90, 0.67, 0.54, 0.43, 0.33, top to 
bottom on the right of (e), bottom to top on the right of (f), and top to bottom on the left of (f). The curves in (g) are for 
r = 0.74,0.60,0.54,0.48,0.43,0.33,0.14 top to bottom on the left and bottom to top on the right. The curves in (h) are for 
different seeds. 



peaks on a log-log scale. The straight line is a linear fit 
for L > 50 giving the estimate 7/1/ = 1.748(11). For the 
estimation oi (3/v we use the values of the order param- 
eter at the estimated critical temperature (Tc = 1.3983). 
As shown in panel (d), on a log- log scale, the linear fit 
provides the estimate fi/v = 0.127(5). Thus, our results 
for the pure 2d BC model at A = 1 are in full agree- 
ment with the findings of Beale [26| and with universal- 
ity arguments that place the pure BC model in the Ising 
universality class. 



IV. PHASE TRANSITIONS OF THE 
RANDOM-BOND 2d BC MODEL 

A. Second-order transition emerging under 
random bonds from the first-order transition of the 
pure model 

Figure [3] illustrates the effects, at A = 1.975, induced 
by bond randomness for different disorder realizations, 
on the dp structure for lattices with linear size L = 30 
(Figs. [3i;a)-(c)) and L = 50 (Fig. M.d)). It is immedi- 
ately seen that the introduction of bond disorder has 
a dramatic influence on the dp structure of the energy 
PDF. The very rough energy PDF of the pure model, 
with the huge oscillations observed in relatively small 
lattices, is radically smoothed by the introduction of 



disorder (Figs. Elja) and (b)) and the energy barrier is 
highly reduced as the disorder strength is increased. This 
barrier reduction effect depends of course on the disor- 
der realization, as can be easily observed by comparing 
Figs.[3ja), (b), and (c), but its main dependence comes 
from the value of the disorder strength r and already for 
r = 17/23 ~ 0.74, the dp structure is completely elim- 
inated. This is clarified by showing the five realizations 
of size L = 50 in Fig. [S^d). Note here that, for this 
value of the disorder strength r = 17/23 ~ 0.74, only 
a very small portion (< 8%) of the realizations at the 
size L = 30 shows a dp structure in the energy PDF, 
but now with a very tiny energy barrier, whereas for the 
same disorder strength, all realizations at L = 50 have a 
single peak energy PDF. We continue our illustrations of 
the disorder effects by showing in Figs, and (f) the 
energy (e) - T behavior and the order-parameter (M) - T 
behavior, for the disorder realization with seed 99. It is 
clear from these figures that, in effect, only the disorder 
strength r = 0.9 resembles a first-order behavior, whereas 
all other disorder strengths resemble the usual second- 
order behavior. In other words, only for very weak disor- 
der strengths the finite-size rounded anomaly resembles 
a discontinuity in the energy and the order parameter. 

From the order-parameter behavior shown in Fig. [3ljf ) 
it should be observed that the low-temperature behav- 
ior, for strong disorder strengths, shows an unexpected 
and rather complex behavior, which is most prominent 
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FIG. 4; Behavior of the random-bond 2d BC model at A = 1.975: (a) Illustration of the clear saturation of the specific heat 
([C']at)) for the random-bond (open symbols) 2d BC model, (b) FSS behavior of five pseudocritical temperatures defined in 
the text for r — 3/5 = 0.6. The lines show a simultaneous fit for the larger lattice sizes L > 50. (c)-(f) Estimation of critical 
exponents v, jv, and jijv for the case of r = 3/5 = 0.6. In all panels the fits have been performed for the larger lattice sizes 
shown (L > 50). 



for r = 1/3 — 0.33 . . .. Figure [3ljg) further clarifies this 
low-temperature effect, by presenting one more disor- 
der realization (seed 1721) for several disorder strengths 
ranging from r = 1/7 ~ 0.14 to r = 17/23 ~ 0.74. It 
is clear from this figure that for these strong disorder 
strengths the ground state of the model deviates appre- 
ciably from the all Si = -1-1 or the all Si — —1 ferromag- 
netic state. Apparently, this deviation strongly depends 
on the disorder strength as shown in Fig. [3jh), where 
the low-temperature behavior of the order parameter is 
presented for three different disorder realizations. From 
this illustration it appears that, in the strong disorder 
regime, the T —^ value of the order parameter averaged 
over the disorder will depend on the disorder strength. 
This observation will have direct relevance to the ferro- 
magnetism enhancement (from quenched bond disorder!) 
and to the conversion of first-order transitions to second- 
order transitions, through the microsegregation mecha- 
nism to be presented and quantified further below. In 
fact, we have fully verified the above observation for a 
small 4x4 square lattice, for which the exact enumera- 
tion of the spin configurations (3^^) is feasible, using 50 



disorder realizations. 

It is evident from the dp structures of Fig. [3] that 
one should avoid working with values of J2 very close 
to J2 = 1 (pure model), since the first-order character- 
istics of the pure model may be very strong and finite- 
size effects will obscure any FSS in relatively small lat- 
tices. We therefore carried out extensive simulations at 
r = 17/23 ~ 0.74 and r = 3/5 = 0.6. Figure g^a) 
contrasts the specific heat results for the pure 2d BC 
model and both disordered cases, r = 17/23 ~ 0.74 and 
r = 3/5 = 0.6. This figure illustrates that the satura- 
tion of the specific heat is very clear in both cases of 
the disorder strength. However, the presented specific 
heat behavior for both disorder strengths is unsuitable 
for any FSS attempt to estimate the exponent ratio a/v, 
as a result of the early saturation of the specific heat. 
However, the early saturation of the specific heat defi- 
nitely signals the conversion of the first-order transition 
to a second-order transition with a negative critical ex- 
ponent a. Furthermore, using our numerical data, we 
attempted to estimate a complete set of critical expo- 
nents for both values of the disorder strength considered 
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here. For r = 17/23 ~ 0.74, our FSS attempts indicated 
that we are still in a crossover regime for the lattice sizes 
studied. On the other hand, for the disorder strength 
r = 3/5 = 0.6, the FSS attempts, using the larger lat- 
tice sizes studied (L = 50 — 100), provided an interesting 
and reliable set of estimates for the critical exponents, 
which seems to satisfy all expected scaling relations. Fig- 
urelU^b) gives the behaviors of five pseudocritical temper- 
atures T[z]^^ corresponding to the peaks of the following 
quantities averaged over the disorder realizations: sus- 
ceptibility, derivative of the absolute order parameter, as 
defined in Eq. ([5]), and first-, second-, and fourth-order 
logarithmic derivatives of the order parameter, as defined 
in Eq. @ . The five lines shown are obtained via a simul- 
taneous fitting of the form T^z]' = Tc + bL~^/^ for the 
larger lattice sizes L > 50. The overall shift behavior is 
very convincing of the accuracy of our numerical method. 
This accuracy is due to the fact that for each realization, 
the WL random walk has been repeated in the first stage 
of the entropic scheme five times, thus reducing signif- 
icantly the statistical errors, which are then further re- 
fined in the second stage of the entropic process. Further- 
more, since these points are derived from the peaks of the 
averaged curves and not from the individual maxima of 
the realizations, they do not suffer from sample-to-sample 
fluctuations and large statistical errors. This good behav- 
ior allows us to estimate, as shown in Fig. [3]Jb), quite ac- 
curately both the critical temperature Tc = 0.626(2) and 
the correlation length exponent ly = 1.30(6). Regarding 
the latter, we shall see below that it agrees with the esti- 
mate obtained via the FSS of the logarithmic derivatives 
of Eq. (|4]) and this will be a very strong indication of the 
self-consistency of our scheme. 

Figures IHc)-(f) give the FSS behavior of the first-, 
second-, and fourth-order logarithmic derivatives of the 
order parameter defined in Eq. ([4]), the magnetic suscep- 
tibility, the absolute order-parameter derivative defined 
in Eq. ([5]), and the order parameter at the critical tem- 
perature estimated in Fig. [3I^b). Figure IDJc) shows a 
simultaneous fit for three moments and for lattice sizes 
L > 50. The resulting value for the exponent ly — 1.35(5) 
indeed agrees with our earlier estimate from the shift be- 
havior in Fig. Bfb) and also fulfils the Chayes et al. in- 
equality V >2/d [331 . Figure mjd) presents the behavior 
of the peaks of the average susceptibility on a log-log 
scale and the solid line shows a linear fit for sizes L > 50. 
The estimated value for the exponent ratio ^/v shown 
in this panel is very close to 1.75 and it is well known 
that this value of the ratio ^/v is obeyed not only in 
the simple Ising model but also in several other cases in 
2d. In particular, it appears that it is very well obeyed in 
the cases of disordered models, including the site-diluted, 
bond-diluted, and random-bond Ising model S 

113, El B E EE 113, H Si, IM Im US, [Ha iHrilf ■ 

Furthermore, it has been shown that is also very well 
obeyed in both the pure and random-bond version of 
the square Ising model with nearest- and next-nearest- 
neighbor competing interactions [13, HE] , as well as the 



case of the second-order phase transition induced by bond 
disorder from the first-order behavior in the g = 8 Potts 
model [HP. FigureUe) is a first estimate for the expo- 
nent ratio P/v — 0.12(1) obtained from the FSS behavior 
of the maxima of the average absolute order-parameter 
derivative [Eq. (O] with the solid line shown being a lin- 
ear fit, again for L > 50. Finally, Fig. [IJf) shows the 
conventional FSS method of estimating the ratio P/v by 
considering the scaling behavior of the average order pa- 
rameter at the estimated critical temperature — 0.626. 
The solid line is a linear fit for L > 50 giving the value 
jS/v = 0.128(5). These latter two estimates are very close 
to the value (3/v = 0.125 and combining the above results 
one finds that the random-bond version of the model ap- 
pears to satisfy the scaling relation 2(3/v + jjv = d. 
Thus, a kind of weak universality appears [47l. [57|. 



B. Second-order transition emerging under 
random bonds from the second-order transition of 
the pure model 

We now present our numerical results for the random- 
bond 2d BC model with A = 1 for disorder strength 
r — 0.6. Bond randomness favoring second-order transi- 
tions, this system is also expected to undergo a second- 
order transition between the ferromagnetic and param- 
agnetic phases and it is reasonable to expect that this 
transition will be in the same universality class as the 
random-bond 2d Ising model. As far as we know, there 
has not been any previous attempt to compare the behav- 
iors of the random-bond BC model and the random-bond 
Ising model. The latter model is a particular case of the 
more general random Ising model (random-site, random- 
bond, and bond-diluted) and has been extensivel y in ves- 
tigated and debatedJTp^^ El, El, El, El,!!, El, 
El, Ei, US [H, m, [lllHirifr^ using renormallzation- 
group and conformal field theories, the marginal irrele- 
vance of randomness at the second-order ferromagnetic- 
paramagnetic transition has been predicted. According 
to these theories, the effect of infinitesimal disorder gives 
rise only to logarithmic corrections as the critical expo- 
nents maintain their 2d Ising values. On the other hand, 
there is not full agreement in the literature for the fi- 
nite disorder regime. Two existing scenarios are mutu- 
ally exclusive: The first view predicts that the specific 
heat slowly diverges with a double-logarithmic depen- 
dence, with a corresponding correlation-length exponent 
V = 1 [H, El, E3, E3l ■ Another scenario predicts a neg- 
ative specific heat exponent a leading to a saturating 
behavior [ll] , with a corresponding correlation length ex- 
ponent V > 2/d. Let us now present the FSS analysis of 
our numerical data. 

Figure E^a) presents again the shift behavior of six 
pseudocritical temperatures, as defined above for the 
pure model, but here using the peaks of the correspond- 
ing quantities averaged over the disorder realizations. 
Fitting our data for the larger lattice sizes {L = 50 — 100) 
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FIG. 5: Behavior of the random-bond (r = 0.6) 2d BC model at A = 1: (a) Simultaneous fitting of six pseudocritical 
temperatures defined in the text for L > 50. (b) FSS of the averaged specific heat peaks. A double-logarithmic fit is applied 
for L > 50. The inset shows a linear fit on a double-logarithmic scale, (c) FSS behavior of the averaged susceptibility peaks on 
a log-log scale, (d) FSS of the averaged order parameter at the estimated critical temperature, also on a log-log scale. Linear 
fits are applied for L > 50. 



to the expected power law behavior, T = Tc + bL^^^", 
we find that the critical temperature is Tc = 1.3812(4) 
and the shift exponent is l/i' — 1.011(22). This lat- 
ter estimate is a first indication that the random-bond 
2d BC at A = 1 has the same value of the correlation- 
length critical exponent as the pure version and there- 
fore as the 2d Ising model. Figure [5Kb) illustrates the 
FSS of the specific heat maxima averaged over disorder, 
[C]*„. Using these data for the larger sizes L > 50, we 
tried to observe the goodness of the fits, assuming a sim- 
ple logarithmic divergence, a double-logarithmic diver- 
gence, or a simple power law. Although there is no ir- 
refutable way of numerically distinguishing between the 
above scenarios, our fitting attempts indicated that the 
double-logarithmic scenario applies better to our data. 
The double-logarithmic fit is shown in the main panel 
and also in the inset of Fig. (S^b). Finally, Figs. [5^0) and 
(d) present our estimations for the magnetic exponent 
ratios j/v and P/v. In panel (c) we show the FSS be- 
havior of the susceptibility peaks on a log-log scale. The 
straight line is a linear fit for L > 50 giving the estimate 
1.749(7) for 7/1/. For the estimation of P/v we have used 
the values of the order parameter at the estimated crit- 
ical temperature Tc — 1.3812. This traditional method, 
shown in panel (d) on a log-log scale, provides now the es- 
timate P/v — 0.126(4). From the above findings, we con- 
clude that, at this finite disorder strength, the random- 
bond 2d BC model with A = 1 belongs to the same 
universality class as the random Ising model, extending 
the theoretical arguments based on the marginal irrele- 
vance of infinitesimal randomness. Most strikingly, it is 
undisputable from our numerical results that the second- 
order phase transitions emerging, under random bonds, 
from the first-order and second-order regimes of the pure 
model, have different critical exponents although they are 
between the same two phases, thereby exhibiting a strong 
violation of universality. We note that, since our bond 



disorder occurs as the variation of the bond stren gths 
that all are in any case non-zero, no Griffiths line [5g| 
divides the paramagnetic phase here. 

Finally, we discuss self-averaging properties along the 
two segments (ex-first-order. Sec. IVA, and still second- 
order. Sec. IVB) of the critical line. A useful finite-size 
measure that characterizes the self-averaging property of 
a system is the relative variance Rx = Vx/[X]^„, where 
Vx = [-'^^lau — [-'^IL: of singular extensive thermo- 
dynamic property X. A system exhibits self-averaging 
when Rx ^ as L 00, or lack of self- averaging (with 
broad probability distributions) when Rx const ^ 
as L — » 00. The FSS scenario of Aharony and Har- 
ris [59! describes self-averaging properties for disordered 
systems and has been validated by Wiseman and Do- 
many [60l . [biI . [s^ I in their study of random Ising and 
Ashkin- Teller models. From these papers, the disordered 
system resulting in Sec. IVB from the introduction of 
bond randomness to the marginal case of the second- 
order transition of the pure 2d BC model is expected 
to exhibit lack of self-averaging. This expectation also 
agrees with our recent study of the self-averaging prop- 
erties of the 2d random-bond Ising model [63|. In the 
current work, the FSS behaviors of the relative variances, 
obtained from the distributions of the magnetic suscep- 
tibility maxima {X = x*)i were observed. Their behav- 
ior clearly indicates that the disordered systems exhibit 
lack of self-averaging along both segments of the criti- 
cal line, since these relative variances show a monotonic 
behavior and are still increasing at the maximum lat- 
tice sizes studied. For A = 1, R^^iL = 90) ~ 0.0011 
and Rx'{L ^ 100) ~ 0.0014, whereas for A = 1.975, 
R^>{L = 90) ~ 0.015 and R^.{L = 100) ~ 0.016. Thus, 
the latter case, i.e., the ex-first-order segment, gives a 
much larger effect, by a factor of ~ 12. A similarly 
stronger lack of self-averaging was observed for the dis- 
ordered system resulting from the case of competing in- 
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teractions on the square lattice, than for the disordered 
system resulting from the marginal case of the simple 
Ising model, in our recent study [s^l- Moreover, the dis- 
cussion in Sec. VIIA in the paper by Fisher [g^I is relevant 
here, explaining the expectation for extremely broad dis- 
tributions near ex-first-order transitions in systems with 
quenched randomness. This discussion points out also 
that ex-first-order transitions may have several v expo- 
nents m, ini and provides the background for 
understanding why our finite-size correlation-length ex- 
ponent obeys the Chayes et al. inequality [soj . 



C. Contrasting random-bond behavior of critical 
temperatures, connectivity spin densities, and 
microsegregation 

In most spin models, the introduction of bond random- 
ness is expected to decrease the phase-transition temper- 
ature and in several cases the critical temperature goes 
to zero at the percolation limit of randomness (r = 
and J2 =0). For less randomness, only a slight decrease 
is expected, if the average bond strength is maintained. 
Indeed, in the second-order regime of the pure 2d BC 
model, A = 1, the introduction of bond randomness 
has slightly decreased the critical temperature, by 1% 
fSec. lIVB]) . On the other hand, in sharp contrast, for the 
same disorder strength r = 0.6 applied to the first-order 
regime of the pure model (Sec. IIV Ap . at A = 1.975, we 
find a considerable increase of the critical temperature, 
by 9%. 

In order to microscopically explain the above observa- 
tion of ferromagnetic order enhanced by quenched disor- 
der, let us define the following connectivity spin densities, 
Qn = {Si)n, where the subscript n denotes the class of 
lattice sites and is the number of the quenched strong 
couplings (Ji) connecting to each site in this class. Fig- 
ure |6| illustrates the temperature behavior of these den- 
sities averaged over 10 disorder realizations for a lattice 
of linear size L — 60. For A = 1.975, it is seen that the 
Si = preferentially occur on the low strong-coupling 
connectivity sites. The Si = ±1 states preferentially oc- 
cur with strong-coupling connectivity, which (1) natu- 
rally leads to a higher transition temperature and (2) 
effectively carries the ordering system to higher non-zero 
spin densities, the domain of second-order phase transi- 
tions. Figure [5] constitutes a microsegregation, due to 
quenched bond randomness, of the Si = ±1 states and 
of the Si = state. We note that microsegregation is 
reached by a continuous evolution within the ferromag- 
netic and paramagnetic phases. A similar mechanism 
has been seen in the low-temperature second-order tran- 
sition between different ordered phases under quenched 
randomness [67| . 

On the other hand, for A = 1 and in the neighbor- 
hood of the critical temperature, the difference between 
the smallest and the largest of the connectivity densities 
is 0.177, whereas for A = 1.975, this difference is 0.449 




0.2 I ' ' ' ' ' ' ' ' ^ 
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FIG. 6: (color online) Temperature behavior of the connec- 
tivity spin densities Q„ defined in the text for A = 1 (upper 
curves) and A = 1.975 (lower curves) for a lattice size L — 60. 
In each group, the curves are for n = 0, 1, 2, 3, 4 from bottom 
to top. 

in the corresponding critical region. Thus, the microseg- 
regation does not occur in the regime of the second-order 
transition of the pure model and the effect of quenched 
disorder is the expected slight (the average bond strength 
is maintained) retrenchment of ferromagnetic order. Mi- 
crosegregation does occur in the first-order regime of the 
pure model where macrosegregation occurs in the ab- 
sence of bond randomness. The result is a local con- 
centration of Si = ±1 states, leading to the enhancement 
of ferromagnetism. The above-mentioned spreads in the 
connectivity-density values are very slowly changing with 
the lattice size. For instance, for L = 20 the correspond- 
ing values are respectively 0.175 and 0.445. 



V. CONCLUSIONS: STRONG VIOLATION OF 
UNIVERSALITY AND WEAK UNIVERSALITY 

In conclusion, the second-order phase transition of the 
2d random-bond Blume-Capel model at A = 1 appears 
to belong to the same universality class as the 2d Ising 
model, having the same values of the critical exponents, 
i.e., the critical exponents of the 2d Ising universality 
class. The effect of the bond disorder on the specific heat 
is well described by the double logarithmic scenario and 
our findings support the marginal irrelevance of quenched 
bond randomness. On the other hand, at A = 1.975, the 
first-order transition of the pure model is converted to 
second order, but in a distinctive universality class with 
ly = 1.30(6) and /3/i/ = 0.128(5). 

These results, on the 2d square lattice, amount to a 
strong violation of universality, since the two second- 
order transitions mentioned in the previous paragraph, 
with different sets of critical exponents, are between the 
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same ferromagnetic and paramagnetic phases. This re- 
sult was also obtained by renormalization-group calcu- 
lations [13] in 2d and 3d that are exact on hierarchical 
lattices and approximate on square and cubic lattices. 
The mechanism in these renormalization-group calcula- 
tions is that the second-order transitions, emerging un- 
der random bonds from the first-order transitions of the 
pure model, have their own distinctive unstable zero- 
temperature fixed point [H, [6^ . 

Furthermore, the latter of these two sets of results sup- 
ports an extensive but weak universality, since the latter 
of the two transitions mentioned above has the same mag- 
netic critical exponent (but a different thermal critical ex- 



ponent) as a wide variety of 2d systems without [H, [T3| 
and with [13, El, HE 113] quenched disorder. 
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